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ABSTRACT 

We introduce a new set of simulations of Milky Way-sized galaxies using the AMR code ART -|- 
hydrodynamics in a ACDM cosmogony. The simulation series is named GARROTXA and follow the 
formation of a halo/galaxy from z = 60 to z = 0. The final virial mass of the system is ^7.4x IO^^Mq. 

Our results are as follows: (a) contrary to many previous studies, the circular velocity curve shows 
no central peak and overall agrees with recent MW observations, (b) Other quantities, such as 
M_*(6xIO^^M q) and R^ (2.56 kpc), fall well inside the observational MW range, (c) We measure the 
disk-to-total ratio kinematically and find that D/T=0.42. (d) The cold gas fraction and star formation 
rate (SFR) at z=0, on the other hand, fall short from the values estimated for the Milky Way. 

As a first scientific exploitation of the simulation series, we study the spatial distribution of the hot 
X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution 
accounting for an important fraction but not all of the missing baryons. An important amount 
of hot gas is also present in filaments. In all our models there is not a massive disk-like hot gas 
distribution dominating the column density. Our analysis of hot gas mock observations reveals that 
the homogeneity assumption leads to an overestimation of the total mass by factors 3 to 5 or to an 
underestimation by factors 0.7 —0.1, depending on the used observational method. Finally, we confirm 
a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems. 

Subject headings: galaxies:formation - methods:numerical - Galaxy: halo 


1. INTRODUCTION 

One of the most challenging problems that hydro- 


works of Evrard 

(1988) 

Hernquist & Katz (19891): Gen 

et al. 

1990(); Navarro & 

Benz 

(|1991); Navarro & White 


hierarchical structure formation scenario, systems that 
look like real disk galaxies; that is, galaxies with ex¬ 
tended disks and present-day star formation rates, disk- 
to-bulge ratios, and gas and baryonic fractions that agree 
with observations. 

It has been long known that disks should for m when 


gas cools and condenses within dark matter halos (White 


& Rees 

1978 

Fall & Efstathiou 

1980 

Mo et al. 

1998| 

conserving angular momentum o 

btained through 

exter- 

nal torques from neighbouring structures I 

Hoyle 

1953 

Peebles|1969|). Yet, the first halo/ 

galaxy hydrodynamical 


20001 inevitably ended up with small and compact disks 
and massive spheroids that dominated the mass. Two 
effects were found to be the responsible of this result. 
First, artificial losses of angular momentum cau sed by 
insufficient resolution and other n umerical effects (|Abadi 


et al.||2003[ Okamoto et al(]|2003| |Governato et al.||2004 


Kaufmann et al. 2007). Second, the dynamical friction 

that transfers the orbital angular momentum of merging 


substructures to the outer halo and causes cold baryons 
(stars plus gas) to sink to the cen ter of the proto-galaxy 


(e.g. Hernquist & Mihos 1995). Once in the center 


baryons that are still in the form of gas are quickly con¬ 
verted to stars. These newborn stars will join the stellar 
component and both will form an old spheroid, leaving 
no g as for the formation o f new disk stars at latter times 
(e.g. Mailer fc DekeH2002 ). Later, the increase in resolu¬ 
tion in the numerical simulations led to several authors 


(e-S- 

Robertson et al.||2004l Okamoto et al.||2005 

Gover- 

nato et al.|2U(J4 

) to succeed in getting more reads 

tic disk 


called zoom-in technique. This technique allows simula¬ 
tors to obtain high resolution by choosing a small region, 
almost always a sphere centered on a specific halo, from 
a rel atively big box, and rerun it with higher resolution 


(e.g. Klypin et al. 20011. To obtain a “MW-like” system 
the authors selected halos with masses similar to the one 

HO^^ M 0 ) and halo s that have 
1.5 


estimated for the MW (^ 
not suffered a major merger since z ■ 


e-g- 


Robertson 


et al. 112004 Governato et al.||2004[ [Ukamoto et al.||2005[ 


Scannapieco et al.|2009||Guecles et al.|2011[|Moster et al.| 

3r et al. 


2014 Vogelsberger et al.||2014p . 

The recent relative success on forming more realistic 
disk galaxies has not been limited to resolution. The 
implementation of new subgrid ph ysics such as an ef- 
ficient supernovae (SNe) feedback (Stinson et al. 2006) 
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and new star formation rec ipes have shown to be im¬ 
portant ( Gnedin et al.|[2009 |. The more effective stellar 
feedback acts by removing; low-angular-momentu m ma¬ 
terial from the central part of (proto)galaxies (e.g. Brook 
et al.|[20T2) while a SF based on the molecular hydrogen 


To solve this missing baryons problem, following t he nat¬ 


ural prediction from structure form ation models d White 
fc Rees 19^ . White&_fte^ 1991^ authors like |Gen 


Ustriker 1 1999); Uave et al. 12001 ); Bregman (2007) pro- 


abundance is more realistic. These improvements make 
simulated galaxies to follow observed scaling relations, 


like t he Tully-Fisher, for the first time (Governato et al. 


2007) and thus relax the tension that existed between 


posed that galactic winds, SiNe feedback or strong AGN 
winds ejected part of the galactic baryons to the dark 
matter halo and circu mgalactic medium (CGM) as hot 


obser vations and AGDM predictions (Governato et al. 


2010). In galaxy formation simulations, gas is accreted 


m hlamentary cold flows that never is shock-heated to the 


2009). 


halo virial temperature (Keres et al. 2005 Dekel et al. 


mation in disks (Keres et al. 12009 Ceverino et al. 

20101). 

h'or instance, authors like Scannapieco et al. (|2009 

: Stin- 

son et al. ( 

2010 I): iPiontek & bteinmetz ( 2 O 1 I); lAgertz 

et al. 

(I 20 II ); Brooks et al. (201l|); h'eldmann et al. 


Gued 

es et al. (2011) obtained realistic rotational 

ly sup- 


gas. Others authors like Mo & Mao (2004) proposed that 
also part of the gas never collapsed into the dark matter 
halos as it was previously heated by SNe of Population 
III. From the point of view of large volume hydrodynam- 
ical galaxy formation simulations, some attempts have 
been made in order to investigate th e presence and de- 
tectability of hot gas halo corona. Toft et ^ (2002) 


presented one of the first studies towards this direction 
however due to the inefficient feedback, the amount of hot 
gas in their simulated halo was too small. Global proper¬ 
ties of hot gas halo cor ona in MW-sized simulations have 
also b een studied (e.g. Guedes et al.|2011 Mollitor et al.j 


MW ones. 

Until recently, simulated disk galaxies with rotation 
curves similar to that of the Milky Way were scarce. 
Most works used to obtain syst ems with too peaked and 


2015). More recently, models that account tor the SINe 
feedback and implement new star formation processes, 
succeeded in obtaining more realistic results whe n com¬ 
pared with observations (Grain et al.||201(3 2013). 


declining rotation curves 

(e.g. Scannapieco et al. 

2012 

Hummels & Bryan 

2012 

). Agertz et al. 

(mini 

rowed 


From an observational point or view the hot gas 
halo corona was proposed to be detecte d by its X-ray 
emissi on. First detections where done by [Forman et al 


that by tine tuning the values ot the Sh' efficiency and (1985) in external galaxies and later on several authors 
the gas density SF threshold parameters it was possible 
to avoid the formation of galaxies with such peaked rota- 


get better detect ions in other systems and studied 


including ours, this is not longer a problem ( 

Aumer et al. 

2013 Mollitor et al. 20151 |Marinacci et a 

i.||20141 |Mu- 

rante et al. 20151 Agertz & Kravtsov 2015| Keller et al. 

2015 Santos-Santos et al. 12016 

, in part due to an added 

“eaHy-teed back” (Stinson et al. 

2013 

) to the thermal SNe 


2008 

lAnderson k5regman||2UlU[ |Bog 

dan ^ Giltanov 

2011 

TTJHi 

. More recently, using Chandra, 
3 and other instruments (e.g. 

XMM-JNewton, 
Hagihara et al. 

2010 

2011 Miller & Bregman 2013a 

it 

was detected a 

large reservoir ot this h 
(e.g. ~6xl0io M 0 in 

ot eas surroundine our galaxy 

Gupta et al. 

2012) and it was 


and/or the inclusion of the AGN feedback. The early 
feedback is attributed to a feedback that is present from 
the very beginning, before the first su pernova explodes . 


and it may include radiation pressure (Krumholz 


2015), 


proposed that it could account for a fraction of the so 
called missing baryons mass. These detections were 
done by the analysis of X-ray absorption lines from 
OVII and OVIII which only exist in environments with 
temperatures between 10® and 10^ K, the MW halo has 


photoheating of the gas due to the ionizi ng radiation of 
massive stars (|Trujillo- Gomez et al.|2015), stellar winds, 
etc. Currently the need of such effects are widely ac- 
knowledge d but their correct implementations are still a 
challe nge (Oman et al.||2015 Gonzalez-Samaniego et al. 
M4|). 

Obtaining systems that resemble the MW open several 
possibilities on the study of galaxy formation and evolu¬ 
tion. For instance, it is interesting to study the origin of 
the galactic halo coronal gas, their properties and how 
they can account for a fraction of the missing baryons. 
The lack of baryons in the universe was reported by 
Fukugita et al. (1998) when they found that cosmolog- 
ical baryon traction inferred from Big Bang nucleosyn¬ 
thesis is much higher than the one obtained by counting 
baryons at redshift z = 0. This problem was confirmed 
when a more accurate value for the cosmological baryon 
fraction was obtained from t he high preci sion data of 
WMAP (jPunkley et a l.j 2009) and Planck (Planck Col- 
laboration et al.|f2014). In tact, the cosmic ratio ilh/ilM 
is 3—10 times larger than the one observed for galaxies. 

^ These authors need to delay or switch off cooling for some time 
once the “regular” SN feedback is switch on to avoid overcooling. 


a temperature of about log T = 6.1—6.4 (Yao & Wang 
2007| Hagihara et al. 2010). The problem of using this 


technique is that these X-ray absorption lines only can 
be observed in the directions of extragalactic luminous 
sources (QSO, AGN, ...) or galactic X-ray emitters 
as X-ray binaries. This observational strategy gives 
a limited information about hot gas distribution and 
position, inside the galactic halo or in the intergalactic 


works (Gupta et al. 

2012 

Miller & Bregman 

2013b 

Gupta et al. 2U14) obtainee 

1 a value tor the total hot 


for about 10—50% of missing baryons. To obtain these 
total hot gas mass from the low number of available 
observations the authors needed to assume a simple 
hot gas density profile. It is important to be aware 
that this assumption c an lead to biased tota l hot gas 
mass values. Recently Faerman et al. (2016) used an 


analytic phenomenological model to study the warm/hot 
gaseous coronae distribution in galaxies. They have 
used the MW hot gas X-ray and UV observations as 
input for their model and found that the warm/hot 
halo coronae may contain a large reservoir of gas. They 
conclude that if metallicity is of about 0.5 solar it can 
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account for an important fraction of the missing baryons. 

Here, we introduce a new set of MW-sized simulations 
with a high number of DM particles and cells (~7xl0®). 
On the other hand, the spatial resolution, the side of 
the cell in the maximum level of refinement, is 109 pc. 
These are highly resolved simulations even for today 
standards. The simulated galaxy ends up with an ex¬ 
tended disk and an unpeaked slowly decreasing rotation 
curve. It is important to mention that the subgrid 
physics used in the present work was shown to produce 
realistic low-mass galaxies, at least as far as the rotation 
curves was concerned (e.g. 

GARROTXA 


Colin et al. 20131. In the 
series of simulations fne temperature 
reached by the gas in the cells, where stellar particles 
are born, is actually higher t han the one obtained in, for 
examp le, the simulations by |Gonzalez-Samaniego et al.| 
(2014) because the SF efficiency factor is higher (0.65 
versus 0.5). This along with the fact that here the gas 
density threshold is lower, 1 cm instead of 6 cm 
makes the stellar feedback in GARROTXA simulations 
stronger. 

Using our simulations we have also studied how the 
hot gas component embedded in the DM halo can 
account for part of the missing baryons in galaxies. 
A novelty of the present paper is the determination, 
in our GARROTXA series of simulations, of the hot 
halo gas distribution in a full sky view which offers us 
an opportunity to detect observational biases in the 
determination of its mass, when only a small number of 
line of sight observations are used. 

This paper is laid out as follows. In Sec. [^we present 
the code and initial conditions we have used. In Sec. [3] we 
introduce the general properties of our MW-sized simu¬ 
lations, whilst in Sec. we present the study of their hot 
gas component. We summarize our conclusions in Sec.[^ 

2. GARROTXA SIMULATION 
2.1. The code 

The numerical simulations we introduce in this work 
(GARROTXA: GAlaxy high Resolution Runs in a 
cOsmological GonteXt using ART) have been computed 



1997 

Kravtsov] 

2003 

). This is the 

0 

t the code introc 

luced 

in IColm et al.l 


(2010) in the context of the formation of low-mass 
galaxies. This vers io n dif fers from the one used by 


Geverino & Klypin (2009), among other things, m 
the feedback recipe; unlike them, we deposit all the 
energy coming from stellar winds and supernovae in the 
gas immediately after the birth of the stellar particle. 
This sudden injection of energy is able to raise the 
temperature of the gas in the cell to several 10^ K, high 
enough to make the cooling time larger than the crossing 
time and thus avoiding most of the overcoolin^ Here 
we have used it to obtain high resolution Milky Way 
sized galaxies. 


The code is based on the adaptive mesh refinement 
technique which allows to increase the resolution selec¬ 
tively in a specified region of interest around a selected 
dark matter halo. The physical processes included in 
this code are the cooling of the gas and its subsequent 
conversion into stars, the thermal stellar feedback, the 
self-consistent advection of metals, and a UV heating 
background source. The used total cooling and heat¬ 
ing rates incorporate Gompton heating/cooling, atomic, 
and molecular hydrogen and metal-line cooling , UV heat¬ 


ing from a co smological background radiation (Haardt & 
Madau|I996), and are tabulated for a temperature range 


of 10^ K < T < 10® K and a grid of densities, metallici- 
ties (from log(Z) = -3.0 to log(Z)=I.0, where Z is in solar 
units) , and redshifts, using the CLOUDY code (Ferland 
I998[ version 96b4). 


et al. 


'I'he star formation (SF) is modeled as taking place 
in the coldest and densest collapsed regions, defined 
by T < Tsp- and rig > nsF, where T and rig are the 
temperature and number density of gas, respectively, 
and Tg^;’ and nsF are the temperature and density SF 
threshold, respectively. A stellar particle of mass m 
is placed in all grid cells where these conditions are 
simultaneously satisfied, and this mass is removed from 
the gas mass in the cell. The particle subsequently 
follows N-body dynamics. No other criteria are imposed. 
In most of simulations presented here, the stellar particle 
mass, m, is calculated by assuming that a given fraction 
(SF local efficiency factor csf) of the cell gas mass, mg, 
is converted into stars; that is, m = csf nig, where esF 
is treated as a free parameter. In MW-sized models 
presented here we have used csf = 0.65, Tg^;’ = 9000 K 
and ngj? = I cm“^, where this last is the density 
threshold in hydrogen atoms per cubic centimeter. As 
shown in Golfn et al. (2010), these values successfully 
reproduce realistic low-mass galaxi es at least as the 
circular velocity was concerned. In Golfn et al. (20I0| 
the authors also found that the reduction ol esF makes 
the temperature reached by the gas in the star forming 
cell to be too low , as a consequence the cooling became 
too efficient. They also show that the strength of the 
stellar feedback recipe used here depends also on the 
value of nsF- the lower its value, the stronger the 
effect of the feedback. We tested that most of the 
overcooling is avoided when ngF is around I cm“® and 
esF = 0.65. Finally, the value Tgi,’ is almost irrelevant 
as long it is set below or close to 10"^ K because this is 
always achieved when the density threshold is reached. 
Incidentally, the Jeans length of an isothermal gas, with 
a density of n//=3 cm“^ and a temperature of 3000 K 
is 638 pc, which is more than four times the length of 
the cell at the highest refinement level at present-day 
(109 pc). This value has been obtained usin g the 
standard perturbative derivation of Jeans criteria ( Kolb ] 
|& Turner) 1990 Shu 1991 Binney & Tremaine 2(JU8|) 
that delivers the expression shown in Fqj^for the Jeans 
length, where Cg is the sound speed and the po the 
mass density of the star forming gas (cold and dense 


^ This temperature is the one that the cell acquires if one as- 
sumes that all the energy injected by the SNe (and stellar winds) 
is in the form of thermal energy. Its exact value will depend on the 
assumed IMF and, filonsp and SNe energy values. 


^ As nsF is 1 cm“^, it is expected that at the moment of the 
stellar particle formation the gas density to be greater than 1 but 
not much greater. Moreover, as these simulations do not have self¬ 
shielding the average gas temperature in the disk lies around the 
few thousands degrees. 
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gas). In Appendix B we show a histogram of the Jeans 
length, in units of the corresponding cell, for the cold gas. 


A,7 = 



( 1 ) 


In this work we have found that these values also 
successfully reproduce a galaxy rotation curve similar 
to the one of the Milky Way at z=0 for galaxies with 
mass similar to the one of the MW. Since stellar particle 
masses are much more massive than the mass of a star, 
typically 10^—10® Mq, once formed, each stellar particle 
is considered as a single stellar population, within which 
the ind ividual stellar m a sses a re distributed according 
to the Miller & Scalo (1979) IMF. Stellar particles 
eject metals and thermal energy through stellar winds 
and Type II and la SNe explosions. Each star more 
massive than 8 Mq is assumed to dump into the ISM, 
instantaneously, 2x10®^ erg of thermal energy; 10^^ erg 
comes from the stellar wind, and the other 10®^ erg 
from the SN explosion. Moreover, the sta r is assumed to 
eject 1 .3 Mq of metals. For the assumed |Miller & Scalo 
(1979) initial mass function, IMF, a stellar particle of 
10^ Mq produces 749 Type II SNe. For a more detailed 
disc ussion of the pr o cesses implemented in the code, 
see |Kravtsov| (|2003|); [Kravtsov et al.| (|2005|). Stellar 
particles dunip energy in the form of heat to the cells 
in which they are born. If sub-grid physics are not 
properly simulated, most of this thermal energy, inside 
the cell, is radiated away. Thus, to allow for outflows, 
it is common to adopt the strategy of turning off the 
cooling during a time to// in the cells w here young 
stella r particles (age < to//) are placed (see Colfn et al. 
2010). This mechanism along with a relatively high 
value of esF allow the gas to expand and move away 
from the star-forming region. As toff can be linked 
to the crossing time in the cell at the finest grid, we 
could see this parameter as depending on resolution in 
the sense that the higher the resolution, the smaller 
its value. In simulations presented here the cooling is 
stopped for 40 Myr. Although it has been demonstrated 
that for spatial resolutions similar to the ones in our 
models this process is u nnecessary when simulating 
dwar f galactic systems (Gonzalez-Samaniego et al. 
2014), a more detailed study is necessary to ensure that 
the same occurs when simulating MW-sized galaxies. 


2.2. Simulation technique and halo selection 

Simulations presented here have been run in a ACDM 
cosmology with fig = 0.3, Ua = 0.7, ^b= 0.045, and h = 
0.7. The CDM pow er spectrum was taken from |Klypin| 
& Holtzman (1997) and it was normalized to cts = 0.8, 
where erg is the rins amplitude of mass fluctuations in 
8 Mpe h“^ spheres. To maximize resolution efficiency, 
we used the zoom-in technique. First, a low-resolution 
cosmological simulation with only dark matter (DM) 
particles was ran, and then regions of interest (DM 
halos) were picked up to be re-simulated with high 
resolution and with the physics of the gas included. The 
low-resolution simulation was run with 128^ particles in¬ 
side a box of 20 Mpe h“^ per side, with the box initially 


covered by a mesh of 128^ cells (zeroth level cells). At 
z=0, we searched for MW-like mass halos (7.0x10^^ Mq 
< Myir < 1.5x10^^ A4 q) that was not contained within 
larger halos (distinct halos). We selected only halos that 
have had not major mergers since z=1.5 and that at z=0 
have not a a similar mass companion inside a sphere 
of 1 Mpe h“^. Other restrictions we have imposed are 
that halos need to be in a hlament or a wall not in 
a void or a knot. After this selection, a Lagrangian 
region of 3 R.uir- was identified at z =60 and re-sample d 
with additional small-scale modes (Klypin et al. 2002). 


The virial radius in the low resolution runs, K,i,ir 5 is 
defined as the radius that encloses a mean density equal 
to Ayir times the mean density of the universe, where 
Ayir is a quantity that depends on Dg, Ha, and z. For 
example, for our cosmology Ayir(z = 0) = 338 and 
Ayir{z = 1) = 203. The number of DM particles in 
the high resolution region depends on the number of 
DM mass species and the mass of the halo. For models 
with four or hve DM mass species this value varies from 
^ 1.5x10® (model G.242) to about 7.0x10® (model 
G.321). The corresponding DM first specie mass per 
particle of modeled galaxies (naoMisp) is given in Tab.[^ 

In ART, the initially uniform grid is refined recur¬ 
sively as the matter distribution evolves. The cell is 
refined when the mass in DM particles exceeds 81.92 (1- 
Fb_[/) mp/toM or when the mass in gas is higher than 
81.92 Ff,// uip/tg, where Ff,// is the universal baryon 
fraction, mp=7.75x10"* Mq h“* is the mass of the 
lightest particle specie in the DM only simulation, and 
fnM and tg are factors that control the refinement 
aggresivitiy in the dark matter and gas components, 
respectively. In our models Fh^(7=0.15 according to the 
cosmology we have imposed and the mean DM and gas 
densities in the box are assumed to be equal to the 
corresponding universal averages. For the simulations 
presented in this work the grid is always unconditionally 
refined to the fourth level, corresponding to an effective 
resolution of 2048® DM particles. A limit also exists for 
the maximum allowed refinement level and it is what will 
give us the spatial size of the finest grid cells. Here we 
allow the refinement to go up to the 11 refinement level 
what corresponds to a spatial size of the finest grid cells 
of 109 comoving pc. Our main model is a re-simulation 
of a MW-like halo of tAyir ^7.4x10** Mq. However 
we also re-simulated other halos with masses in between 
13.9x10** Mq and 2.0x10** Mq. We have used these 
simulations to hnd how the general properties of the 
galactic system change with mass. For each one of these 
models we have also made simulations changing initial 
parameters such as csf, "^sf^ resolution or number of 
DM mass species. This set of models helps us to assess 
how the hnal results depend on the initial conditions. A 
more detailed study of the effects of changing the sim¬ 
ulation para meters can be find in Gonzalez-Samaniego] 
et al. (2014). Tab. fl] summarizes the properties of the 
main MW-sized simulations presented here. 

Although our MW-sized simulations are realistic, some 
physical processes are not well implemented in the code 
we have used. For instance, the AGN feedback, the 
molecular cooling down to 200 K, the formation of molec- 
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ular H 2 clouds, a deep chemical treatment, the radiation 
pressure or the photoionization from massive stars. How¬ 
ever, the last two are now implemented in a new version 
of the code and a new set of simulations is being gener¬ 
ated. 

In this work we study MW-sized systems, however it is 
not our goal to reproduce all observed properties. Here 
we define a MW-sized system as the one that has a total 
mass, Vmax and disk with similar mass and scale lengths 


Robertson et al.||2004[ |Governato et ai.||2004[ [Okamoto 

et al.||2005[ |Scannapieco et al.||20091 |Guedes et al. 

2011 

Moster et al. 120141 

Aumer et al.|2013||Stinson et al. 

2013 

Vogelsberger et al. 

2014p we have selected a DM halo that 


in the high resolution run at z=0 has no other massive ha¬ 
los inside a 800 kpc sphere and that have had a quiescent 
recent assembling history. We note that our initial condi¬ 
tions do not reflect the environment in which the MW is 
embedded but this is certainly not needed if we just aim 
at producing a system with a non-negligible disc com¬ 
ponent. This latter is almost always achieved under the 


in last 10 Gyr or so (e.g. 

Moster et al. 120141 

Aumer et al. 

2013 

Stinson et al.|2013 

Vogelsberger et al.| 

i'hese 


tails of each realization will depend on selected subgrid 
physical parameters, environment and assembling histo¬ 
ries. 


3. GARROTXA MW-SIZED RUNS 
3.1. General properties 

In this work we focus mostly in three MW-sized 
models (G.3 series) that differ only in the aggresivity 
on the refinement criteria. Parameters that control 
the refinement aggresivity (Jdm^ fg) take the following 
values: (l.,l.) for the less aggressive model, (4.,8.) for 
the intermediate and (8., 32.) for the more aggressive. 
These models are the ones we have found to be in 
a better agreement with observed properties of the 
MW. The spatial resolution of all our models is high, 
around 109 pc. Model G.321 is the one with a lightest 
refinement, G.322 has intermediate conditions and 
G.323 is the one with hardest refinement conditions 
(harder means that the code uses a lower mass threshold 
for the rehnement). All parameters we present in this 
section are summarized and compared with recent 
MW-sized models and with its observational values in 
Tab. [2 In the table we see that there are parameters 
which are very sensitive to a change in the refinement 
setting as, for example, the SFR(z=0) that changes 
from 0.02 to 0.32, G.322 and G.323, respectively. Yet, 
most parameters do converge as can be seen in Tab. 0 
in particular, circular velocities profiles agree within 5% 
(see Appendix A for more information). 

To ensure numerical convergence we also generated a 
low resolution model. Comparison between models can 
be found in the Appendix A. 


3.1.1. Morphology 

At z~0 our three G.3 models are massive spiral galaxies 
with several non-axisymmetric structures in the disk such 
as bars or spirals (see Fig.f^andj^. Their assembling his¬ 
tory is quiet after z=1.5 and theJust major merger occurs 



[kpc] [kpc] [kpc] 


Fig. 1.— Face on (first column) and edge on view (last two 
columns) of stars in our simulated Milky Way sized galaxy. Each 
row correspond to a different stellar age. From top to the bottom, 
11-13.467 Gyr, 10-11 Gyr, 9-10 Gyr, 4-9 Gyr and 0-4 Gyr. 
Color code indicates stellar age. 


at z=3. In Fig. we show the face on (first column) and 
edge on (last two columns) views of stars of our model 
G.323 at z=0. For the other two models the main picture 
is similar but with a smaller number of stellar particles. 
Each row correspond to a different stellar age (see figure 
caption). In this figure it can be observed that young 
stellar component (0—4 Gyr) is distributed in a flat disk 
structure. On the other hand, older stellar populations 
(4—10 Gyr) are also distributed into a disk structure but 
in these cases not as flat as the younger. We argue that 
disk scale height and length depends on the age of the 
stellar population and that the former is higher and the 
latter is smaller when the po pulati on gets older, this is¬ 
sue will be addressed in Sec. 3.4.2| Another interesting 
property that can be easily observed in the edge on views 
of our models is that stellar population has an evident 
flare that is more evident for older populations. The face 
on view of the younger stellar population shows the pres¬ 
ence of rings, spiral arms and also a young bar, in the 
disk. Spirals, rings and the young bar can be better ob¬ 
served in models G.322 and G.323 where the number of 
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Fig. 2.— Edge on (top) and face on (bottom) views of a young 
stellar population (0—7 Gyr) of model G.323 at z=0. Both Panels 
span 50 kpc in the x-axis. The y-axis span 22 kpc in the top 
panel and 50 kpc in the bottom. Color code indicates stellar age, 
qualitatively. 

stellar particles is higher (see Fig.[^. The young bar has 
grown from disk particles that have ages between 0 and 
8—9 Gyr. The old stellar population (11—13.476 Gyr) is 
distributed in a bar/ellipsoid component and it has been 
originated in the last major mereer at z=3 (Valenzuela 
et al. in preparation). In Fig. ^we show how the gas 
component as function of temperature (see figure cap¬ 
tion). In this figure is easy to see that gas distribution 
is not isotropic. Gas is distributed in different regions 
of the system depending on its temperature: cold gas is 
present in the young stellar disk region and hot gas fills 
the out-of-plane region and is embedded in the DM halo. 
It is also interesting to see how cold/intermediate gas 
has a warped structure that is also marginally observed 
in the stellar component. 


3.1.2. Stellar, gas and dark matter virial mass 

Here we define the virial radius {vyir) as the one where 
the sphere of radius Vyiy encloses a mean density 97 times 
denser than the critical density {pcru) of a spatially flat 
Universe /9crit=3H^(z)/(87rG). We have used value 97 as 
it is the value derived from the spherical top-h at col¬ 
lapse model for A GDM at z=0 for our cosmology (Bryan 
& Norman 1998). This virial radius definition is bor¬ 
rowed from structure growth theory and then its use 
is not appropriated when on e wants to define a phys¬ 
ically meaningful halo edge (Cuesta et al. 2008 Zemp 
2014|). To avoid this problem we also give the prop¬ 


erties for another commonly used definition which is 
i' 200 j the radius that encloses a mean density equal to 
200 times Pcrit- Using the first definition we have ob¬ 
tained that the virial radius at z=0 is lyiy ~230 kpc in 
our three main models and that the mass enclosed in 
this radius is 20—7.61x10^^ Mq. This value for 

the Myiy falls well inside the observational ran g e that 
is M„ir=0.6— 2.4x 10^^ Mq by Xue et al. (2008), Kafle 



Fig. 3.— Face on (left) and edge on (right) views of gas density 
at z=0 in model G.323, as function of its temperature. From top 
to the bottom, 10®-10^ K, 3xl0®-10® K and 10^-3x10® K. In 
color we show the gas temperature. 


et al. (2012), Boylan-Kolchin et al.| (2013) and Kafle 


et al. (2014). I'he total mass enclosed at r2oo=175.6 kpc 


is M 2 oo= 6 . 71 — 6.90x 10^^ M©. Total virial mass is dis¬ 
tributed in dark, stellar and gaseous matter as follows: 
M£,m=5.86-6.79x10“ Mq, M*=6.1-6.2x101° Mq and 
MgQs=l-73—2.70x 10^° Mq. The baryonic fraction of 
our G.3 models is Fb[/=0.104—0.121. These values 
are 19—31% smaller than the universal value for the 
adopted cosmology which is Fb_j/=0.15. Inside lyir, and 
for the model G.321, we have a total number of par¬ 
ticles of N(otai=7.52x10°, where N£iMisp=7.12x 10° and 
N»=3.94x 10° and a total number of gas cells of 2.Ox 10°. 
For models with a more aggressive refinement the number 
of stellar particles and gas cells grow up to 2.34x 10° and 
6.8x10°, respectively. All dark matter particles inside 
virial radius belong to less massive DM specie (Isp) and 


have a mass of 9.25x10^ Mq. Star particles have masses 
between ~10° and 1.2x10° rboca 

are summarized in Tab. [T] 


3.1.3. Circular velocity curves 
Circular velocity curves of our main simulations, com¬ 
puted using the ^GM^y/r approach, are shown in Fig.|^ 
top left panel, as a solid black line. Circular velocity 
curves have also been computed using the real galactic 
potential by using Tipsy package. Results obtained us- 
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ing both techniques are in good agreement. Fig. also 
shows circular velocity c urves from other stat e of the art 
MW-si zed simulations : Klyp in et al. ((2002) B 1 model 
(blue) , | Mollitor et al. (2015) model (magenta), [Guedes 


et al. (| 2011 ) FKib simulation (green). As a compari¬ 


son w e show in red the analytical fit to thelSofue et al. 


MW rotation curve data, presented in Pichardo 


et al.| (2003). Solid lines are the total circular velocity 
curves while the DM, stellar and gas components are 
shown as long dashed, dotted and short dashed lines, re¬ 
spectively. In the top right panel we show, first, G.321 
velocity rotation curves of stars and gas, computed as the 
mean of the tangential velocity in rings centered to the 
galactic center (cyan solid and dashed lines), and the to¬ 
tal circular velocity curve obtained using Nc=\fGM<^ 
approach. As expected, gas follows the circular velocity 
while stars are affected by the asymmetric drift. In this 
panel we also compare stellar rotation curve of our sim¬ 
ulation (cyan solid line) with values from MW observa¬ 
tions. We show rotation velocity curve estima tions from 


blue h orizontal-branch halo stars in t he SPSS (Xue et al . 
2008) (cyan and magenta dots), from Sotue et al. (|2009[) 


(red d ots) and from observations of Lopez-Corredoira 
(2014) (green dots). In the bottom panel we compare 


th e recent dat a for t he circular velocity curve obtained 
by (Reid et al.| (|2014|), using high mass star forming re- 
gions, with the total circular velocity curve of our G.321 
model. As can be seen in the figure our G.3 models 
are in a a very good agreement with the most recent 
observational data. We also show the recent hypothet¬ 
ical gaseous_j;otationciu've of the Milky Way obtained 


{ 2015|) after correcting for bar non- 
t IS important to mention that we 


by Chemin et al. 
circular motions, 
only show results for G.321 model because circular ve¬ 
locity curves of models G.322 and G.323 do not differ 
significantly. The peak of z=0 circular velocity of our 
models is reached at Rpeafe ~5.69 kpc with a value of 
Vc(Rpeafe)=237.5— 243.8 km s“^, the value at a standard 
solar radius (R 0=8 kpc) is Vco=233.3—239.8 km s“^. 
Also the ratio between circular velocity at 2.2 times disk 
scale radius (V 2 . 2 ) and circular velocity at r 2 oo (V 200 ) 
of our simulations (V 2 , 2 /V 2 oo "^1-9) is inside the obser¬ 
vational range for the MW that is 1.671 q'24 ~ 20 

( Xue et ^(2008 Dutton et al.||2010 ). 

An interesting exercise has been to compare circ ular ve¬ 
locity curve s of our simulations with the one of |Guedes| 
et al. (2011) model. We have easily seen that our simula- 


tions have a bit higher velocity c urves and do not pre sent 
a peak in the inner regions as in Guedes et al. (2011). In 
simulations this internal peak is usually associated with 
the presence of a massive bulge in the central parts of the 
simulated galactic system. On the other hand, a similar 
peak i n the Vn is shown in observational work of |Sofue] 


et al. (2009 ). It is in discussion whereas this peak de¬ 


tected in Sofue et al. (2009) observations is a signature 
of non-circular motions inside the bar or of a massive 


bulge ( Duval fc Athanassoula|1983 Chemin et al. ]M5t . 
Non-circular motions need to be computed to solve such 
question. We plan to undertake this work in a near fu¬ 
ture following the one started by Chemin et al. (2015). 
Finally, despite of the rotation and circular velocity 
curves we introduce here fall well inside observational 
ranges it is important to take into account that nowa¬ 


days observational uncertainties are still high. 













































































^vir 

M.^ 27 ' 

M* 

Mpas 

^hotgas 

^^warmgas 

^^coldgas 


[kpc] 

[Mq] 

[M©] 

[Mq] 

[M©] 

[M©] 

[M©] 

G.321 

230.1 

7.33x1011 

6.1x101® 

2.70xl0i® 

1.22xl0i® 

5.66x10® 

9.34x10® 

G.322 

230.1 

7.20x1011 

6.1x101° 

1.73xl0i® 

0.98xl0i® 

5.98x10® 

1.52x10® 

G.323 

230.1 

7.61x1011 

6.2xl0i° 

1.96xl0i® 

1.32xl0i® 

4.54x10® 

1.86x10® 

ERIS 

239.0 

7.90x1011 

3.9xl0i° 

5.69xl0i® 

3.60x101® 

14.20x10® 

6.70x10® 

MollitorB 

234.0 

7.10x1011 

5.6x101® 

7.96x101® 

— 

— 

— 

MW obs. 

- 

l.OibO.30x1012 

4.9-5.5x101° 

- 

- 

- 

7.3-9.5x10® 


r200 

M 200 

Fb.u 

raoMlsp 

I^*min 




[kpc] 

[Mq] 

[M©] 

[M©] 

[M©] 

[M©] 

G.321 

175.6 

6.84x1011 

0.120 

0.93x10® 

0.24x10'! 

120.0x10! 


G.322 

175.6 

6.71x1011 

0.109 

0.93x10® 

0.12x10'! 

120.0x10! 

— 

G.323 

175.6 

6.90x1011 

0.107 

0.93x10® 

0.12xl0! 

120.0x10! 

— 

ERIS 

175.0 

6.60x1011 

0.120 

0.98x10® 

0.63x10! 

0.63x10! 

2.0x10! 

MollitorB 

176.5 

6.28x1011 

0.191 

2.30x10® 

4.50xl0! 

4.50xl0! 

— 

MW obs. 

- 

- 

- 

- 

- 

- 

- 


Npart 

NoAflsp 

N. 

Npas 

Resolution 

CPU time 

CODE 


[10®] 

[10®] 

[10®] 

[10®] 

[pc] 

[IQ! h] 


G.321 

7.52 

7.12 

0.39 

2.0 

109(lcell) 

2.5 

ART + hydro 

G.322 

7.95 

7.08 

0.86 

5.0 

109(lcell) 

4.3 

ART + hydro 

G.323 

9.44 

7.10 

2.34 

6.80 

109(lcell) 

9.2 

ART + hydro 

ERIS 

15.60 (-I-SPH) 

7.00 

8.60 

3.0 

120(e) 

160.0 

GASOLINE (SPH) 

MollitorB 

6.06 

3.91 

2.15 

2.50 

150(lcell) 

— 

RAMSES (hydro) 

MW obs. 

- 

- 

~100000.0 

- 

- 

- 

- 


c=r„ir/rs 

Rd 

,young 

^z,old 

^^hotgas / Mi;i7’ 

ax 

SFR (z=0) 



[kpc] 

[pc] 

[pc] 



[M© yr-i] 

G.321 

28.5 

2.56 (4.89/2.21) 

277 (exp) 

1356 

0.016 

-0.62 

0.27 

G.322 

26.8 

3.20 (5.26/2.26) 

295 (exp) 

960 

0.013 

-0.83 

0.02 

G.323 

26.9 

3.03 (3.19/3.31) 

393 (exp) 

1048 

0.017 

-0.67 

0.32 

ERIS 

22.0 

2.50 

490 (sech^) 

- 

0.017 

-1.13 

1.10 

MollitorB 

56.5 

3.39 

— 

— 

0.046 

-4.54 

4.75 

MW obs. 

21.Ul« 

to 

b 

1 

300±60 

600-1100±60 

- 

- 

0.68-1.45 


Vc0(R=8 kpc) 

I^pea/c 

Vc(Rpeafc) 

V2.2/V200 

nsF 

Tsf 

CSF 


[km s“l] 

[kpc] 

[km s“l] 


[cm“®] 

[10® K] 


G.321 

239.8 

5.69 

243.8 

1.90 

1.0 

9 

0.65 

G.322 

233.6 

5.69 

237.5 

1.85 

1.0 

9 

0.65 

G.323 

233.3 

5.69 

237.8 

1.93 

1.0 

9 

0.65 

ERIS 

206.0 

1.34 

238.0 

1.66 

5.0 

30 

0.10 

MollitorB 

233.0 

9.63 

233.0 

— 

2.7 

3 

0.01 

MW obs. 

221±18 

- 

- 

1 1 1 

-^•-^-^-0.20 

- 

- 

- 


Box 

Zini 

no 

S2a 


Ho 

<^8 


[Mpc h“^] 





[km s ! Mpc !] 


G.321 

20 

60 

0.30 

0.70 

0.045 

70 

0.80 

G.322 

20 

60 

0.30 

0.70 

0.045 

70 

0.80 

G.323 

20 

60 

0.30 

0.70 

0.045 

70 

0.80 

ERIS 

90 

90 

0.24 

0.76 

0.042 

73 

0.76 

MollitorB 

20 

50 

0.276 

0.724 

0.045 

70.3 

— 

MW obs. 

- 

- 

- 

- 

- 

- 

- 


TABLE 1 Parameters of GARROTXA321 (G.321), GARROTXA322 (G.322) and GARROTXA323 (G.323) simulations at z=0. is 
assumed to be equal to rgy. Mass related parameters: is total mass inside Vyir. M* and ^-re total stellar and gas mass inside 

Tyir. Hot, warm and cold gas mass ? ^warmgas) '^coidgas) ^-re gas mass inside Tyir with temperature T > 3x10^ K, 3x10^ K > 

T > 3x10^ K and T < 3x10^ K, respectively. If is the baryonic fraction m£)jvfisp is the mass of a single particle 

belonging to the first DM mass specie, "i^*Tnin/max is the lowest and highest mass of stellar particles and mspH is the mass of SPH 
particles. Structural parameters: c is the DM halo concentration parameter defined as ryir/vs, Rd is the disk scale length, ^z,young/old is 
the disk scale height of young (0—0.5 Gyr) and old (4.0—11.0 Gyr) stellar populations and ax is the hot gas density power-law exponent 
{photgasM ^ r“^). Kinematical parameters: Vc© is the circular velocity at solar position (R© = 8 kpc), Rpeafc is the radius in where 
circular velocity curve reaches its highest value (Vc(Rpeafc))- Initial condition parameters: ngp, Tsp and are the star formation 
density and temperature thresholds and the star formation efficiency. Box is the simulated box size in Mpc h~^, Zini is the initial redshift 
and ^bi H p and crs define the assumed cosmol ogical model. Compar ison with other state of the art simulations is also shown 

(ERIS simulation IlGuedes et al. | 20l']l and model B inIMollitor et al.|(|2 0 1^L Obse r ved v a lues for the MW have been obtained fro m 
iFerrierel(|20011l: iFlynn et al.| d^UUbll: [jJu~ et al.|||2006 ulHain mer et al.l JijuOyII: iXue et alT](|2008ll: lJuric et al.N2008ll: iDutton et al.|||2010|| and 

[Katie CT al. l j2U14| |, see text for more details on reterences. 


3.2. Dark matter component 


As we have mentioned in the last section all dark mat¬ 
ter particles inside virial radius are particles of the first 
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- G.321 (this work) 

- Pichardo et al. 2003 

• 

Lopez-Corredoira et al. 2014 

- Klypin etal. 2002 (Bl) 

^rotyStoT 


^SDSSi 

- Guedes et al. 2011 

- - - , y 

' Tot,gaa 

1 

^SDSS, 

- Monitor etal. 2014 (B) 

1 Sofue et al. 2009 



-^^-'- 300| -^^^- -n - 



R[kpc] R[kpc] 



^total 


(G.321) 

— Chemin et al 2015 
i I Reid et al. 2014 


OL 


10 

R[kpc] 


15 


20 


Fig. 4.— Top-left: Circular ve locity c urve of our simulated Milky Way sized galaxy model G.3 21 fblackUGuedes et al.|||2011| (green), 
model B1 in |Klypin et al.| < |2002| | (bl ue), |Mollitor et al.|||2015[) (magenta) and analytical model in [Pichardo eraZ] ( |2UU3l l fred). 'the figure 
shows the contribution to the circular velocity Vc=\/GM<r/r of dark matter (short dashed curve), stars (long dashed), gas (dot-dashed) 
and total (solid curve) mass components. Top-right: Total circular velocity curves of the same models shown at left, observational data 
points and gas and stars rotation curves of our model G.321 (cyan dashed and solid lines, respectively). Data point s come from tw o 
realizations of the rotation curv e of the Milky Way from observations of blue horizontal-b r anch halo stars in the SDSS (jXue et al.|2008|l, 
in cyan and magenta dots, from |L6pez-Corredoira| (|2014||, in green, and fro miSofue et an( 2009ll, in red. Bottom: Observational circular 
velocity curve obtained using massive young star terming regions of the MW l|Keid et al.|2UM I. I he total Vc of our model G.321 at z=0 is 
shown as a black solid line. In green we show the hypothetical gaseous rotation curve ot tUe Milky Way obtained by |Chemm_et_^ 
after c orrecting the non- circular motions in th e velocity profile of the Milky Way inferred with the tangent-point method by jKalberla^f^L] 
< |2005| |. In grey we show [Chemin et al.| | [20T^ curve la range. 
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Fig. 5.— Average dark matter (red dots) and hot gas (T > 3x 10^ 
K, blue dots) density profiles at z=0. The upper solid yellow line 
show the best NFW profile fit to the dark matter mass distribu¬ 
tion, in the range from 5 kpc, outside bar-bulge region, to r^ir. 
The lower solid yellow line show the best fit power-law profile to 
the hot gas mass distribution between 5 and 20 kpc. The DM best 
fit NFW profile is characterized by a large halo concentration pa¬ 
rameter c=28.5 as the dark matter halo contracts in response to 
the condensation of baryons in its center (adiabatic contraction). 
In the inner region it can be clearly observed that no DM core is 
present. The hot gas best fit power-law profile gives a slope of 
qx=- 0 . 62 . 


DM specie, the less massive one. In Fig. [^we show the 
dark matter density profile as a function of radius (red 
circles) and the best fit of NFW profile for the model 
G.321 (upper yellow solid line). We avoid the central re¬ 
gion (R < 5 kpc) in the fit of the NFW profile to avoid 
perturbations derived from the presence of baryons (adia¬ 
batic contraction). The best NFW fit gives a large value 
for the halo concentration that is c=r„ir/rs=28.5, 26.8 
and 26.9 for G.321, G.322 and G.323, respectively. These 
values fall i nside the observational range presented in 
Kafle et al. (2014) that is 21.1^g'*3®. No core is observed 
m the central region o f our G.3 models. T he halo spin 


parameter as defined in Bullock et al. (2001) is A’=0.019, 
0.017 and 0. 022. These A’ values a lso fall near the range 
obtained by Bullock et al. (2001) after analyzi ng more 


than 500 sin iulated halos (see hgure 1 and 2 in Bullock 


et al. (2001)). 


3.3. Gas component 

In our models G.321, G.322 and G.323 the total gas 
mass inside the virial radius is Mgas=2.7, 1.73 and 1.96 
xl0^° M©. In G.321, 9.34x10® Mq of gas mass is in the 
cold gas phase (T < 3x10^ K). In G.322 and G.323 cold 
gas mass is 1.52 and 1.86x 10® M©. The cold gas mass in 
G.321 is comparable to the total mass of the molecular, 
atomic and warm ionized medium inferr ed for the MW 
that is ~7.3—9.5 X10® Mq (Ferriere|2001 ). For the rest of 
models presented here the amount of cold gas is smaller 
than the observed value. This z=0 cold gas mass reduc¬ 
tion is consequence of an increase on the star formation 
at early ages, that in its turn is provoked by the higher 


aggresivity in the refinement criteria. This increase on 
the early star formation is the only significant effect we 
have observed when changing the refinement criteria and 
it is consequence of resolving a larger number of small 
dense regions in where the SF criteria is accomplished. 
At z=0, the SF criteria is satisfied only in the inner disk 
regions (R < 7 kpc). This spatially limited star forma¬ 
tion accords with the young stellar distribution shown in 
Fig. Hand also with the decrease of the Star Formation 
Rate (SFR) at z=0 that can be seen in Fig. 10 


The hot gas mass (T > 3x10® K) is around 1.22, 0.98 
and 1.32x10®° M© for models G.321, G.322 and G.323, 
respectively. The definition of hot gas we use in this 
work is observationally motivated: total hot gas m ass has 
to be infer red from ionized oxygen observations ([Gupta 


et al.|20i2 ) and this is only possible when gas has opacity 


in the X-ray region (T > 3x10® K). In Fig. I we show 
the gas distribution as function of temperature in model 
G.323. Each panel show the distribution of gas at differ¬ 
ent temperatures (see figure caption). As can be seen in 
this figure cold gas is placed in disk region while warm- 
hot is embedded within the dark matter halo. Gontrary 
to the standard assumption, hot gas do not follow the 
dark matter radial distribution. This can be seen in 
Fig. [5] where we show the dark matter density distri¬ 
bution (red) and the hot gas density distribution (blue). 
The hot gas density distribution inside r = 100 kpc can 
be fitted by a power law (piiotgas (r) oc r“^). After fit¬ 
ting data from models G.321, G.322 and G.323 we have 
obtained that hot gas density distribution scale factors 
(ax) of these models are -0.62, -0.83 and -0.67, respec¬ 
tively. These results are slightly higher than values used 
to fit analytical models of hot gas in dark matter halos (- 
0.9 in Anderson & Bregman|2010 ) and the ones obtained 
in oth er similar models like ERla (-1.1 in Guedes et al.j 
2011|) . The hot gas density distribution power law ht tor 


model G.321 is shown in Fig. (bottom yellow straight 
line). 

In Sec. [^ we present our first results on the study of 
hot gas distribution and a comparison with observational 
values. 


3.4. Stellar component 

The total stellar mass in the virial radius of our 
G.3 models is M*=6.1—6.2x 10®° M©, that is com¬ 
parable with values estimated fo r the Mil ky Wav 
(4.5— 7.2x10®° M fn in Flynn et al. (2006) and [Licquia 
& Newman (2015)). However, if we observe tte relation 
between stellar mass and total mass (see Fig.[^ derived 
from our models (blue dots) and also from observations 
of the MW (shadowed region) we see that most of them 
fall quite above the M*/M„ir predicted by cosmological 
theories. The same result is observed when using data 
from o ther recent MW-size d simulations li ke ERIS simu¬ 
lation (Guedes et al.||2011) or the ones of Mollitor et al. 
(2015). In this work we have not studied this mismatch 
deeply, we leave it for the future. 


3.4.1. Stellar disk and stellar spheroid 

Simulated galactic systems in our main models (G.3 
series) have both stellar disks and stellar spheroidal com¬ 
ponents (halo and/or bulge). To study properties of stel¬ 
lar galactic disks we need to know which stars belong to 
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it. To undertake the selection process we have used an 
extension of the k i nemat ic d ecomposition propo s ed by 
Scannapieco et al. (2009). In Scannapieco et al. (2009) 
the authors used the ratio between the real stellar par¬ 
ticle angular momentum perpendicular to the disk plane 
(jz) and the one obtained assuming stellar particles fol¬ 
low a circular orbit (computed from the circular velocity 
curve, jc)- The circular velocity curve has been com¬ 
puted using both, GM/r approach and the real galac¬ 
tic potential; we have found that results do not depend 
on the Vc computation technique. Using jzjjc ratio it 
is possible to distinguish particles that are rotationally 
supported {jz=jc for disk stars) from ones that are not 
{jz < jc for spheroid stars). Here we have decided to 
make additional cuts to improve this technique. We have 
done complementary cuts in metallicity, in the vertical 
coordinate \z\, in the angle between the rotation axis of 
each particle and the vector perpendicular to the plane 
(cos(q;)) and in the distance to the galactic center (R). 
To avoid kinematic biases we have made cuts that do 
not involve kinematics first. We have started the selec¬ 
tion process by making the metallicity cut; next, we have 
made a cut in the vertical coordinate, and later in radius, 
in cos (a); finally we have used jzj jc condition. Some of 
these restrictions we have imposed require from a previ¬ 
ous knowledge of the disk plane position. To find this 
position we have u sed an iterative ge ometrical approach 
(pages 9 to 11 in Atanasijevic 1971). In this approach 
we end up with the plane that minimizes the cumulative 
distance of all particles to it. In the first iteration of this 
process we only use stars that are inside a 20 kpc sphere 
centered within the center of mass of the lightest DM 
mass specie particles. After this first iteration we only 
use stars that fall nearby the newly defined plane (i.e. 
\znew\ < 5 kpc). We have checked that after few itera¬ 
tions we get a plane that coincides with the disk plane 
defined by young stars and cold gas (see bottom panels 
of Fig. Fig. [ 2 ] and Fig. [^ . 

After analyzing the stellar distribution of our G.3 mod¬ 
els in both, Cartesian and phase spaces, we have found 
conditions that better distinguish disk from spheroid 
particles. In our main models these conditions are: 
jz/jc > 0.55, l^l < 3.5 kpc, log(ZIa/ZIa 0 ) > -0.5, where 
Zla© = 0. 00178 and it is the ir on solar abundance ac¬ 
cording to|Asplund et al.| (|2009|), cos(q!) > 0.7 and R < 
25 kpc. All particles inside vvir that do not accomplish 
one or more of these restrictions are considered spheroid 
particles (halo and/or bulge particles). To account for 
thick disk particles we relax the log(ZIa/ZIa 0 ) and jz/jc 
conditions; that is, stellar particles will also belong to 
the disk component if they fulfill all previous conditions 
except that they have log(ZIa/ZIa 0 ) < -0.5 and jz/jc 
deviates by less than 0.75 from jz/jc peal^ The re¬ 
sult of this selection can be seen in Figj^where we show 
the stellar mass fraction as a function w jz/jc for model 
G.321, the result is similar for the other models. In this 
figure we show how well we trace the two main stellar 
components of the model: one centered at jz/jc=^^ i-e. 
rotationally supported, and another at jz/jc=Q- 

The mass of the rotationally supported component 
(stellar disk) is around = 1.82, 2.17, 2.21xlO^°M 0 

^ For example, if the distribution peaks at 1.0 then a particle 
with jz/jc = 0.3 would belong to the disk. 


and that of the spheroidal one (bulge -I- halo) is 
= 4.29, 3.93, 3.79xl0i° M 0 for models G.321, 
G.322 and G.323, respectively. Comparing these re¬ 
sults with observations we find that the mass of the 
spheroidal component in our models (bulge + halo stars) 
is high er than the upper limit obtained by [Flynn et al. 
(2006) which is Mb < 1.3x10^° Mf?; or than th e most 


recent observations presented by Valenti et al. (2016) 
(Mb=2.0±0.3x 10^° M 0 ). The spheroids of our galactic 
systems are ^ 2 times heavier than the one of the MW. 
On the other hand, disks in our G.3 models ar e less 
massiv e thaii the one estimated for the MW by [Bovyj 
& Rix (2013) (^4.6x10^° M 0 ). As a consequence, the 


kinematic disk-to-total ratio (D/T) decomposition, idiak. 


computed as M*disfe/M* (Scannapieco et aL||2010 2015) 


where M*disfc is the mass of stellar particles with j^/j 
> 0.5, results in a value (~0.42) that is significantly 


Scannapieco et al. 2011 

However, as pointed 

by Scannapieco et al. 

12010 

2015 

) the use of a 


kinematical decomposition results in a lower D/T than 
those obtained using a photometric decomposition, as 
the one determined in MW observations. For instance, 


Scannapieco et al. (2010) the authors obtained a set 


of MW-sized galaxies with a kinematical idisk ~ 0.2 and 
they demonstrated that f disk becomes 0.4—0.7 when 
using a photometric decomposition. 

Is the presence of this massive spheroidal component 
due to a large stellar concentration in the central region 
of our simulated galaxies, a common result in earlier hy- 
drodynamical zoom-in simulations? As can be seen in 
the rotation curves shown in Fig. fflthis is not the case 
(the curves do not present a centrm peak). We have an¬ 
alyzed this spheroidal component and we have detected 
that it is mainly a triaxial structure that w as form ed in 
the major merger occurred at z = 3 (see Sec. 3.1.1|). This 


triaxial structure has a low density in the disk region and 
thus it is not a classical massive bulge. 

3.4.2. Disk properties 

With the stellar disk/spheroid selection process we 
have found that disk component has 2.0, 6.0 and 
18.7x10^ particles in each one of our main G.3 models. 
We have also computed the mean volume and surface 
density of the disk, as function of radius. We have 
found that a single or a double exponential power law 
can be fitted to the data depending on the age of the 
selected population. In Fig. we show the total disk 
surface density (black) as function of radius of the 
model G.321 at z=0. We have fitted two exponential 
power laws (red and blue) to the surface density curve. 
Results from these fits are that scale length in the inner 
regions (2—6 kpc) is Rd=4.89 kpc and in the outer ones 
(6—12 kpc) Rd=2.21 kpc. If we fit a single exponential 
to the whole radial range from 2 to 12 kpc we obtain that 
Rd=2.56 kpc, result that is in agreeme nt with values 
obtain ed for the MW (e.g. 2. 3T0.6 kpc in H amme r et al. 
(|2007|) or 2.15±0.14 kpc in |Bovy & Rixj (|2U13|) using 
G-type dwarfs from SEGUI^. It is important to note 
that, as can be seen in Fig. a small concentration of 
stars is present in the centralregions (R < 2 kpc). This 
concentration is caused by the presence of a young bar. 
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Fig. 6.— Comparison between stellar/baryonic fraction 
theoretical predictions and the ones obtained 
from observations and simulations. Blue items show M* 
values in our models G.321 (circles), G.322 (triangles) and G.323 
(squares), using total stellar m ass (lar ge items) and only disk stel- 


disk selection groces^ Ubserved ran^ 

e for the M* in the 

Milkv Wav Galaxv lIFlvnn et al. 120061 

Xue et al.|2008| |Licquia &:| 

|Newman|2015| |Gibbons et al.|2U14[| is 

region. Red dot snow tne M*/M,,,v 

shown as a grey shadowed 

value obtained in model B 


the ERIS sirnulation ||GuMes et al.||2011|l. |Behroozi et al.| ||2013|l 
observational M»/M„iV curve interred trom observations at z=U.i, 
with no distinction between blue or red galaxies, and its errors, is 
shown as dashed black line . As a red dashed line, results ob tained 
for the relation in|Gonzalez-Samaniego et al.|(|2014|| when 

using a semi-empirical model, me blu e dasned line c orresp onds 
to the Ml,/Muir relation calculated by |Baldry et al.|||2008l. As 
can be seen it becomes a challenge for ACUIVI to ht its theoretical 
predictions with observations. 


In Fig. we show the volume density as function of the 
distance to the plane (z) of three disk stellar populations 
split by age. We have selected a young (0—0.5 Gyr), 
intermediate (0.5—4.0 Gyr) and old (4.0—11.0 Gyr) 
populations. The density have been computed using 
vertical bins of Az=0.2 kpc and only stellar particles 
that are in between R=2 and R=9 kpc. We have found 
that the scale height (hj.) of each population, when 
fitting a decreasing exponential to the vertical density 
distribution in the range 0.1 kpc < \z\ < 2.0 kpc, 
is h2=277 pc for the young, hz=959 pc for the in¬ 
termediate and hj;=1356 pc for the old populations. 
These values have been obtained from G.321 model, 
results from models G.322 and G.323 can be found in 
Tab. [2 These scale height values are compatible with 
ones expected for the MW: observations show that the 
vertical distribution of stars in the MW can be fitted by 
two exponential w ith a scale heights of h,:=3 00±60 pc 


( Juri c et al. 20081 and h^=600—1100±60 pc (Du et al. 


20Ub|. Uur values also follow the relation proposed by 


Yoachim fc Dalcanton| (2006) that argue that disk scale 


height in edge on galaxies increases in disks following 
zo «2h^=610(Vc/100) km s"!. 



Fig. 7. — Stellar mass fraction as function of the orbital cir¬ 
cularity parameter jz/jc, describing the degree of rotational sup¬ 
port of a given stellar particle (black: all particles; red: spheroid; 
blue: disk). The stars in a centrifugally supported thin disk man¬ 
ifests itself in a sharply peaked distribution about unity. We have 
plotted only stellar particles inside disk region {\z\ < 0.5 kpc and 
r < 15 kpc) of model G.321 at redshift 0. 



Fig. 8. — Stellar disk surface density profile as function of cylin¬ 
drical radius R (black solid line). The solid blue and red lines 
show the best-fit exponential profiles to the inner and outer disk 
components, respectively. The inner component (2—6 kpc) has an 
R^=4.89 kpc and the outer (6—12 kpc) an Kd = 2.21 kpc. If we 
fit the whole disk density profile from 2 to 12 kpc we obtain an 
Kd = 2.56 kpc. All results have been obtained analyzing model 
G.321 at redshift 0. As a gray s hadowed reg i on we show MW 
surface density profile obtained by |Bovy &: Rix( | |2013| |. 


3.4.3. Star formation history (SFH) 


In our simulations stars are being formed with a 
SFR^=o=0.27 Mq yr-\ 0.18 Mq yr"! and 0.41 Mq yr-\ 
for models G.321, G.322 and G.323, respectively. These 


values are s maller than ones inferred by [Robitaille] 


& Whitney (2010) using Spitzer data, for the MW 


(SFR=0.6 8—1.45 Mq yr ^) or the one from Licquia & 
Newman (2015) that is 1.65±0.19 Mq yr“^. In Fig. 10 




















































































13 



z [kpc] 

Fig. 9.— Vertical density profile of young (blue), intermediate 
(green) and old (red) disk stars in the model G.321 at redshift 0. 


top panel, we show the SFR as function of the redshift 
for the total (black), spheroidal (red) and disk (blue) 
components. We have computed the SFH of each stellar 
component using the formation time of their correspond¬ 
ing stellar particles (saved by the code) inside ryjr at z= 0. 
We have followed the strategy described in Sec. 3.4.1 to 
distinguish between disk and spheroid stellar particles. 
We also show, as a shadowed region, the predicted SFH 
for MW like halos derived from semi-analytical models 
that combine stellar mass functions with merger histo¬ 
ries of halos |Behroozi et al. 2013). The peak in the 
SFH of our G.3 models occurs at slightly early ages than 
the one predicted for a MW-sized galaxy. In the bottom 
panel we plot the total stellar mass as function of red¬ 
shift. Some important results that can be appreciated 
from this figure are, first, that the spheroidal compo¬ 
nent is being build at high redshifts while disk starts its 
formation at around z=2.5 (~11 Gyr from the present 
time), just after the last major merger that was at z=3. 
Second, it is also important to note that the star for¬ 
mation of spheroidal component decrease quickly after 
z^2.4 and become negligible at around z=0.5. The disk 
SFR also decreases fast in the last time instants except 
for some short periods of star formation that are a conse¬ 
quence of the accretion of small gaseous satellites. This 
reduction in the star formation when reaching z=0 is a 
consequence of the reduction of cold gas available for star 
formation. Several processes can lead to a reduction of 
the cold gas mass fraction. First, cold gas can be drasti¬ 
cally consumed at higher redshifts due to a non-realistic 
implement ation of physical pa rameters that controls star 
formation (Liang et al. 2016). When such problem ex¬ 
ists a large amount ot old stars are present at z=0 in 
the inner disk region. In our G.3 models, as can be seen 
in Fig. an old spheroid is present in the disk region, 
however this old stellar population is not big enough to 
account for the total reduction of the observed cold gas 
component. Second, if stellar and SNe feedback is too 
efficient, an important amount of cold gas can be heated 
up and then to become unavailable for star formation. 
In this second scenario what we would expect is that 
the fraction of hot gas increase with time, a behavior 
that is not observed in our G.3 models. Finally a de- 


Sphenjid 

Disk 

Total 

Guedes et al. 2011 
Monitor et al. 2014 
MW predicted values 
MW SFR(z=0) 



Fig. 10.— Star formation history (SFH), top, and total mass 
history, bottom, of all star particles identified within virial radius 
at z=0 in model G.321. Black, blue and red filled dots show the 
total, disk and spheroid star formation rates (top panel) and stellar 
mass (bottom panel), respectively, as a function of redshift. Grey 
shadowed region ind icates the nredicted SFH for a MW-like halo 
with M97=10^^ Mq (Behroozi et al.|2013|l. G reen and red big dot s 
show total SFR at z= U ot models presented by|Guedes et al.|||2011|l ; 
|Mollitor et al.|||2015|l (B), respectively. The blue tuied re gion snows 
the z=U MW bFK range derived from observations by |Robitaille| 
|fc Whitneyl | |2010l l. 

crease in the cold gas inflow due to inhomogeneities in 
the circumgalactic medium (GGM) can also be a cause 
of the reduction of the SFR. In this last hypothesis we 
would observe only a small amount of cold gas present 
around our system, or falling from filaments, in the last 
Gyr. We argue that in our models the reduction of the 
amount of cold gas is a consequence of a combined effect 
of the first and the last hypothesis as we observe that 
an old spheroid is present and also that there is not a 
large amount of cold gas falling from filaments at z=0, 
however we leave the deeper study of such processes for 
the future. 

In ART, the metallicity Z is divided according to which 
kind of SNe produce the metals: ZII or Zla, if they are 
produced by SNH or SNIa, respectively. The metallicity 
Z is then Z = Zla+ZII. Alpha elements, as oxygen, are 
mostly produced by SNH while iron is mostly produced 
by SNIa, we thus can assume that the abundance of the 
former are proportional to ZII while the iron abundance 
is proportional to Zla. 

3.5. Are our G.3 models MW-like systems? 

Although it is not the aim of this paper to obtain 
a realistic model of the MW galaxy, nevertheless the 
models we present here (i.e. G.321, G.322 and G.323) 
can be considered at least MW-sized as they reproduce 
several of the MW observed properties. On the other 
hand, it is also evident that some parameters of these 
models fall far from observations of the MW. Below, 
we give an overview of the parameters that match MW 
observations and of those ones that do not. 
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3.5.1. MW-like properties 

Our study is focused on a galactic system, the G.3 se¬ 
ries, that forms spiral arms and a galactic bar. Its assem¬ 
bly history is quiet after z=1.5 and there are no major 
mergers afte r z=3, something that is a lso expected for 
our Galaxy (Forero-Romero et al. 2011). Its total mass 
is inside the observ a tional range proposed for the MW 
by Xue et al. (|2008 1 ; Boylan-Kolchin et al. (20091; Kafle 
et al. (2012 2014) that is 1.0 ± 0.4 x 10^^ M©. The stel¬ 
lar mass inside virial radius ( 6.1 — 6.2 x 10 ^° M©) falls 
also inside the observational range for the MW that is 


also inside tne obser vational range tor tne MW tnat is 
4.5 — 7.2 X 10^° Mq ( Flynn et al.|2006 [Licquia fc New¬ 
man Finally, the cold gas mass iiiside virial radius 

in model G.321 (9.34 x 10® ) also ma tch the observed 

values (7.3 — 9.5 x 10® Mq in Ferriere (20011). Several 
structure parameters of our G.3 models are also within 
the MW observations. This is the case for the dark mat¬ 
ter halo concentration parameter (c) that in our models 
takes values from 26.8 to 28.5 while obs ervations predict 
it should be 21 .llg^ 3 ® ( Kafle et al.|2014 ). Hot gas power- 


law profile index (ax), defined in Sec. 
our G.3 models is in between -0.62 and -tTK 


3.3 


and that in 
is also close 


to value proposed by Anderson & Bregman (2010) which 
is -0.9. Disk scale lengths and heights of our galactic sys¬ 
tems (G.3 series) also coincide or are close to observed 
values. The former in our model s is 2.56—3.2 kpc and 
observed values are 2.3T0.6 kpc in Hammer et al. (20071 


and 2.15±0.14 in Bovy & Rix (2013). Ubserved disk scale 
height of young and old sta rs is 300 ± 60 pc in|Juric et al.| 


(2008) and 600—1100 pc in Du et al. (2006), respectively, 
and in our models we have obtained for young and old 
populations, 277—393 pc and 960—1356 pc, depending 
on the model. Finally, also several kinematical parame¬ 
ters agree with the MW observations; the rotation curve 
of our simulated galaxies, for instance, roughly match 
observations (see Fig. [^. Some examples of such agree¬ 
ment are the r atio V 2 . 2 /V 200 th at is 1.9 in our models 
and 1.67tQ;24 


Xue et al. 


(2008), and the Vc(R 0 ) that is 


233.3—239.8 km s ^ in this work and 221±18 in |Kopos'w 


et al. (2010) or 236±11 in Bovy et al. (2009). 


3.5.2. Parameters that do not match observations 

Among the parameters that do not agree with the 
observed values of the MW is the SFR and the bulge 
to disk ratio. The cold gas fraction in the G.323 and 
G.323 models also do not match observations, they ac¬ 
count for only 1.52 and 1.86 x 1 0® M^ w h ile ob served 


Ferriere (2001). On 


values are 7.3 — 9.5 x 10® Mq 
the other hand, the SFR in our models ranges from 
0.18 to 0.41 Mq yr“^ values belo w the observational 
rang e 0.68 — 1.45 Mq yr“^ from Robitaille & Whit- 
(|2010|). Finally, we have presented in Sec. 1314.1 
the stellar spheroid-disk decomposition and shown that 
a massive spheroid exists in our models. Such massive 
structure has not been observed in the MW. The to¬ 
tal mass of the spheroid (bulge-|-halo) and the stellar 
disk in our G.3 models are 3.79—4.29 xlO^® Mq and 
1.82—2.21 xlO^® Mq, respectively. Moreover, observa¬ 
tions show that the bulge and the disk mass of our Galaxy 
are <2.3x10^® Mq (Flynn et al. 2006 Valenti et al. 2016) 


and ^4.6x10^® Mq (Bovy & Rix 2013), respectively; 
that is, the spheroid of our models is at least two times 
more massive than the bulge of the MW. 

4. MISSING BARYONS PROBLEM AND THE X-RAY 
LUMINOUS HOT GAS 

4.1. The missing baryons problem 

It was long been known that the cosmological baryon 
fraction inferred from Big Bang nucleosynthesis is much 
higher than the one obtained by cou nting baryons at 
redshift z=0 (e.g. Fukugita et al. 1998). However, it was 
recently, with WMAF and Flarick hig h precision data 
of the cosmic microwave background (Dunkley et al.j 
'2009 Planck Gollaboration et aL]|2014 ), that the cosmic 
baryonic traction was well constrained and thus the 
lack of baryonic mass in galaxies became evident. This 
data revealed that the cosmic ratio between baryonic 
and total matter (Db/Dxf) is 3—10 times la rger than 
the one obs erved in galaxies. For instance, [Hoekstr^ 
et al. (2005) found a baryonic fraction for isola ted spiral 
galaxi^ of 0 .056 and for elliptical of 0.023 while [Dunkleyl 


et al. (2009) presented a cosmic baryonic fraction from 
WIVIAP 5-year data of about 0.171±0.009. 

Studies like Cappi et al. (2013); Georgakakis et al. (2013) 
tried to solve this missing baryon problem proposing 
that galactic winds, SNe feedback or strong AGN winds 
ejected baryons to the circumgalactic medium (GGM). 
Others proposed that most of the gas never collapsed 
into the dark matter halos as it was previously heated 
by SNe of Population HI (e.g. Mo fc M^|2004 |, this is 
known as the preheating scenario. In this last scenario 
missing baryons are still in the extragalactic warm-hot 
intergalactic medium (WHIM), following filaments. 
Thus, it has become clear that studying the GGM is 
necessary not only for finding the missing baryons, 
but also to understand its properties in the context 
of galaxy formation and evolution models. To further 
constrain galaxy evolution and find out which one of the 
proposed theories solves the missing baryon problem it 
is needed to measure the amount of warm-hot gas phase 
present in the GGM and also its metallicity. Results 
of such measures give information about the heating 
mechanism: hot gas with low metallicities is indicative 
of a Pop HI SNe preheating while a more metallic gas 
suggest that it has been enriched in the disk via SNe 
and stellar winds. 

Supporting the first hypothesis, the idea of very strong 
feedback driving galactic-scale winds, causing metal- 
enriched hot gas to be expelled out of the star-forming 
disk possibly to distances comparable to or beyond the 
virial radius was generally accepted. This process also 
lies in the heart of obtaining realistic simulations of dis k 
galaxies (Guedes et al. 2011 Marinacci et al. 2014). 
By this mechanism gas is temporally placed tar trom 
star forming regions delaying the star formation and 
preventing the formation of an old bulge in the center of 
the system. 

From the analysis of OVH and OVHI X-Ray absorption 
lines observed towards the line of sight of extragalactic 
sources, several authors performed estimations of the hot 
gas mass in the MW halo, suggesting that an important 
part o f the missing baryons is located in this hot gas 
phase dGupta et al.||2012[ |2014|). They concluded that 
the warm-hot phase of the GGM is extended over a large 
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region around the MW, with a radius a mass around 
100 kpc and 10^° M©. However, these estimations 
are limited to a few directions towards extragalactic 
sources (bright QSOs) and as a consequence issues like 
the homogeneity or isotropy of the hot gas distribution 
possibly affecting the mass estimation are hard to be 
addres sed by them. Other studies like [Feldmann et al.| 
(2013) also suggest that the presence of this hot gas 
phase can explain part of the isotropic gamma-ray back¬ 
ground observed in Fermi Gamma-ray Space Telescope. 


4.2. Hot gas component in G.321 
4.2.1. Spatial distribution 

The amount of hot gas in our G.3 models is in between 
M/iot = 0.98—1.32x10^° Mq. This hot gas that has a 
temperature above 3.0x 10° K (has opacity in the X-rays) 
is embedded in the dark matter halo but mostly outside 
the stellar disk. Aiming to perform a fair comparison 
with observations we estimated the hot gas Hydrogen 
equivalent column density and emission measure (EM)[^ 
as: 


Ni{ = / njjdl 

J los 

EM = f nidi 

J los 


( 2 ) 


( 3 ) 


where hh is the hydrogen particles density, dl is an 
element of the path in the line of sight and Ug is the 
electron density. 

In Fig ^ we show a full sky view of hot gas column 
density distribution in our model G.321. This figure has 
been obtained computing the hot gas column density 
from a position that is at 8 kpc from the galactic center, 
inside the simulated galactic disk, resembling the Sun 
position, and assuming arbitrary azimuthal angles. We 
have obtained column density values that fall near the 
observational ra nges, if a solar met allicity is assumed 
(see Tab. |3| and Gupta et al.||2012|). However, when 
assuming lower nietallicities, observations are consistent 
with values that are above the ones measured in our 
simulations. This change can give us information about 
the absorption of the local WIMP. We also conclude that 
the distribution of hot gas inferred from all estimations 
is far from homogeneous. In order to quantify the hot 
gas anisotropy we followed two strategies. First we 
computed the amplitude of the first spherical harmonics 
(Y™ from l,m=0 to 5 ) and later we also com puted 
the filling factor as defined by |Berkhuijsen| (|1998|). Our 
results show that the dominant spherical harmonic is 
Y°, that is an indicator of the dipolar component, with 
a small contribution of several high order components. 
The asymmetry of the hot gas distribution with respect 
to the disk plane suggest some degree of interaction with 
the extragalactic medium. The computation of the mean 
filling factor from line of sights distributed all along the 
sky gives as a value of f„ ~0.33±0.I5. This value for the 

^ In this work we are not showing the dispersion measure values 
(DM) as they only differ by a fa ctor of two from the N^f ones. N// 
values are presented in Fig. |13| Although not presented h ere, DM 
values in ou r models are of the same order as the ones in[Guedes| 
|et al.||20lT||. 


filling factor also suggests that hot gas is concentrated 
in a few regions of the sky. It is also impor tant to 
mention that in our G.3 models (see Fig. and 13) we 
do not detect a hot gas thick disk component. This is 
important as it is under discussion what is the contri¬ 
bution of the hot ga s thick disk to the observed X-ray 
emission/absorption (Savage et al. 2003). Additionally 
to the spatial distribution we study the metallicity and 
velocity components of hot gas cells. Metallicity and 
velocity fields provide us information about how the hot 
gas component interacts with the extragalactic medium. 
In Fig. we show the metallicity distribution of the 
hot gas as seen from the Galactic center. It is evident 
that gas metallicity distri but ion does not depend on the 
azimuthal angle (0) (Fig. 14 top panel) while it is clear 


the dependency on the vertical latitud (</>) (Fig. 14 
bottom panel). It is also noticeable a low metallicity 
component located at the northern hemisphere that is 
absent at the southern, where metallicity is in aver age 
much higher. It is important to note that in Fig. 14 


bottom panel, we can distinguish a clear break in the 
hot gas distribution at low latitudes. This hole coincides 
with the position of the stellar/cold gas disk. 

In order to shed some light on the bimodal origin in the 
metals distribution, we have analyzed a set of three color 
maps showing the projected metallicity and v eloc ities 
along each one of the principal planes (see Fig. 15 and 
Fig. 16). Fig. 15 (X-Y plane, i.e. disk plane) shows that 
in general hot gas in the dis k p lane rotates following the 
disk rotation direction. Fig.[^(X-Z plane, top and Y-Z, 
bottom) shows a more coinplex scenario with several 
vertical motions and metallicity gradients, it becomes 
clear that both low and high metallicity vertical flows 
are present in the hot gas halo. It is also clear from 
this figure that a bimodality in metallicity exists in the 
vertical direction and also that low-metallic hot gas flow 
has a mean motion that brings it from outside to inside 
the system, i.e. it is falling from the IGM. High-metallic 
hot gas flows observed in the southern hemisphere have 
a slower and more irregular motion than ones in the 
northern. While it is easy to understand and interpret 
why low-metallic hot gas is falling inside the system, 
and why in some cases hot metallic gas departs from 
the stellar disk region (SN feedback, stellar winds...), it 
is not so obvious why some clouds of such hot metallic 
gas, with sizes of tens of kpc, behaves differently. We 
suspect that this last high-metallic component could 
be associated with a small gaseous satellite that is 
passing through the system at z=0 (see Fig 17). To 


confirm all these hypothesis we have made t he s ame 
velo city -metallicity maps as the ones in Fig. 15 and 
Fig. for several snapshots back to z = 0.5. From the 
analysis of these snapshots we have confirmed that in 
some instants hot gas with enhanced metallicity departs 
from the disk due to feedback of stellar component. 
We have also observed that low metallicity gas appears 
coming from the CGM. Finally, at around z=0, we have 
seen that a gaseous satellite approaches to the main 
galaxy perturbing hot gas flows. The perturbation of 
the hot gas distribution by accretion of small satellites 
is an interesting process and it will deserve a more 
detailed study. The case of LMG has been used as a 
constrain to the hot gas distribution (Salem et al.|2015 ). 
Finally we have also compared the disk plane gas mean 
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Fig. 11.— Total hot gas mass inside = r200 as function of 
total virial mass. Black dots: Simulations from GARROTXA run, 
their properties are presented in Tab. Blue dots: Simulations 
outside GARROTXA run, parameters from these simulations are 
not presented in Tab. 


metallicity with that one of the hot gas in the halo and 
we have found that both components have similar mean 
metallicities that is around -0.64 dex (except for the 
high and low metallicity flows discussed above). 

A conclusion from the metallicity and kinematic analysis 
is the presence of hot gas flows from the IGM, either 
from the low metallicity environment of from sinking 
metal enriched satellites. We also conclude that the 
disk has a low hot gas concentration as a consequence 
of SNe explosions and stellar winds pushing h ot gas up 
to the halo. Such results are coherent with the IDekel fc] 


Birnboim (2006) proposal of cold gas flows feeding halos 
with Muir <10^ M(7). 


4.2.2. Halo virial mass and the total hot gas in galaxies 

In Tab. and Fig. E] we present results suggesting a 
possible correlation between halo virial mass and the to¬ 
tal hot gas in galaxies. After analyzing MW-sized models 
but also few more massive and some lighter models, we 
have found that Myir correlates with total hot gas mass. 
Blue dots are simulated galaxies that come from a work 
in preparation. These simulations will be discussed in an 
extended study in preparation, but we show them to il¬ 
lustrate the correlation persistence regardless of the sub¬ 
grid physics. Their corresponding host halos were drawn 
from a 50 Mpc h“^ boxside and chosen to be relatively 
free of massive companions inside a sphere of 1 Mpc h~^. 
A similar result was presented in Crain et al. (2010). In 
their work they found instead a correlation between X- 
ray luminosity and virial mass. More specifically they 
reported two linear correlations, one for M^^^ < 5x10^^ 
and another for Myir > 5x10^^. Our study is mainly 
focused on MW-size galaxy systems i.e. their low halo 
mass fit, consist ently we have found a simple linear cor¬ 
relation (see Fig. 11). Although absolute values obtained 
i n our work can not be directly compared with the ones 
ir Crain et al. (2010) as we are using Mhotgas and they 
are using Lx, it hnally conhrmed this linear correlation 
will provide a new constrain to virial mass in galaxies, 
including the Milky Way. 


4.2.3. Total hot gas mass: are results from observations 

biased? 


As previously mentioned it is challenging to infer 
the total hot gas mass in the MW by using the avail¬ 
able observational data. Matter in the CCM is hot 
(T = 10^-10^ K) and tenuous (n ~ lO-^-lQ-^ cm-^). 
Plasma in these conditions couples with radiation mainly 
through electronic transitions of elements (C, N, O) in 
their He-like and H-like ionization stages. The strongest 
of these transitions falling in the soft X-rays. Due to 
the extreme low densities and small optical length of 
the CCM, detection of the absorption/emission features 
produced by this material has proved particularly 
challenging with the still limited sensibility of the 
gratings on board Chandra and XMM-Newton. As a 
consequence, observations in only a few number of line 
of sight directions have been useful to observe X-ray 
absorption in quasar spectra. In this scenario it is clear 
that simulations can play an important role in the study 
of hot gas as they let us to explore the complexities in a 
simulated galactic systems. Including such information 
enable us to study how observational techniques can 
lead to biased results for the total hot gas in the galactic 
system. 

With such aim we have located our mock observer 
at 8 kpc from the galactic center (i.e. similar to 
solar position) in the galactic plane and at a random 
azimuthal angle, inside our MW-sized simulations. For 



like oxygen we have assumed that we are observing all 
hot gas, and present the column densities measured in 
each line of sight. As it is common in the literature (e.g. 


Cupta et al. 

2012 |Bregman & Lloyd-Davies 

2007 

), we 

have assumed a spherical uniform gas distri 

bution tor 


the hot gas halo. As a first attempt, we have estimated, 
from the column density values at each direction, the 
total galactic hot gas mass, assuming a path length of 
L = 239 kpc = lyir. We have found that using the 
spherical uniform approach, the total hot gas mass is 
systematically overestimated (see Tab. last column). 
This is a consequence of that the density distribution in 
our model is n ot u niform but it decreases as a power-law 
(see Fig. and [T^. However, this calculation is not very 
useful, as it requires an a-priori knowledge of the CCM 
path length (L) which is observationally unknown. 


A popular technique among observational studie s (i.e. 
Cupta et al. 2012[ [legman fc Lloyd-Davies||2007 ) is to 


take advantage ol the two observables for the hot gas 
in the CCM: the column density (Njj) and the emission 
measure (EM). This technique uses these two quantities 
to derive the total hot gas mass with no need of imposing 
an a-priori mean density or optical path length. We show 
results from using this technique in Tab. both for the 
optical path length, the mean density and the total hot 
gas mass. In this case, as can be seen in the table, re¬ 
sults show an underestimation of the total hot gas mass 
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Model 

Halo # 

z 

M 200 
[loll Mol 

r 200 

[kpc] 

Mstar 

[10“ Mq] 

[lOio Mol 

^^hotgas 

1109 Mol 

^^coldgas 

[lO® Mol 

<^SF 

DMsp 

^re/ 

G.240 

1 

1.0 

8.34 

183.05 

1.31 

3.64 

23.4 

12.6 

0.10 

3 

11 

G.240 

1 

0.25 

11.3 

203.09 

1.63 

4.78 

30.2 

17.2 

0.10 

3 

11 

G.240 

1 

0.0 

13.0 

211.71 

1.75 

6.86 

35.2 

18.8 

0.10 

3 

11 

G.241 

1 

0.67 

8.94 

186.89 

1.40 

3.60 

21.7 

11.3 

0.50 

3 

11 

G.242 

1 

1.0 

7.65 

179.28 

1.04 

2.76 

15.0 

10.6 

0.70 

4 

10 

G.242 

1 

0.25 

10.2 

198.91 

1.12 

3.63 

28.1 

6.33 

0.70 

4 

10 

G.242 

1 

0.0 

11.4 

207.35 

1.15 

5.23 

30.8 

1.63 

0.70 

4 

10 

G.250 

3 

1.0 

2.55 

128.57 

0.14 

0.95 

1.52 

6.24 

0.70 

4 

11 

G.250 

3 

0.25 

4.37 

151.82 

0.18 

2.63 

5.88 

16.0 

0.70 

4 

11 

G.250 

3 

0.0 

5.30 

164.98 

0.20 

2.23 

6.39 

12.4 

0.70 

4 

11 

G.260 

4 

1.0 

0.982 

92.20 

0.03 

0.75 

0.287 

6.17 

0.5 

3 

11 

G.260 

4 

0.25 

1.49 

106.64 

0.04 

0.73 

0.476 

5.79 

0.5 

3 

11 

G.260 

4 

0.0 

1.76 

113.50 

0.04 

0.80 

0.572 

5.49 

0.5 

3 

11 

G.320 

5 

1.0 

4.20 

145.64 

0.60 

8.83 

3.86 

4.10 

0.60 

4 

11 

G.320 

5 

0.25 

6.01 

168.45 

0.63 

1.18 

9.41 

1.26 

0.60 

4 

11 

G.320 

5 

0.0 

6.67 

175.59 

0.64 

1.46 

10.0 

0.405 

0.60 

4 

11 

G.321 

5 

1.5 

3.70 

139.71 

0.54 

1.33 

4.98 

6.17 

0.65 

5 

11 

G.321 

5 

0.67 

5.42 

161.59 

0.58 

2.66 

8.84 

8.99 

0.65 

5 

11 

G.321 

5 

0.25 

6.37 

171.98 

0.61 

1.78 

8.04 

7.91 

0.65 

5 

11 

G.321 

5 

0.0 

6.84 

175.59 

0.61 

2.17 

10.4 

9.34 

0.65 

5 

11 


TABLE 2 

M 200 vs. HOT GAS MASS INSIDE R 2 OO IN OUR SET OF N-BODY PLUS HYDRODYNAMICS SIMULATIONS. We ALSO SHOW THE MOST RELEVANT 
CHANGES IN THE INITIAL PARAMETERS, FROM THE ONES USED IN THE MODEL G.321 (SEE TAB.[^F0R THE DEFINITION OF THE PARAMETERS). 



r [kpc] 


Fig. 12.— Spherically averaged hot gas density (T > 3x10® K) 
as function of radius for the model G.321 at z = 0. 

by factors ^0.7—0.1, result that is independent on the 
observed direction. Again this result is a consequence 
of that the real density distribution of our model is not 
uniform but decreasing with radius as a power law. 

As a conclusion we state that using both, column density 
plus imposed optical depth (independently of using the 
value from a single observation or the mean) or a combi¬ 
nation of column density and emission measure, we are 
not able to obtain the real hot halo gas mass, when as¬ 
suming it is isotropically distributed in the galactic halo. 
Nevertheless, these measurements, particularly the ob¬ 
servational method in the literature using both Nh and 
EM are useful to get an order of magnitude estimate of 
the total hot gas mass of the halo. The work is in progress 
to find the density profile that will allow us to get the 
best hot halo gas mass estimation, from observations. 


4.2.4. Total hot gas mass: accounting for the missing 
baryons 

The baryonic fraction in our simulations correspond¬ 
ing only for stellar and cold gas mass, falls in between 
0.084 and 0.096. These numbers ar e far from the cosmic 
baryon fraction Fb^c/=0.171±0.009 ( Dunkley et al.|[2009 



1 

[deg] 

b 

[deg] 

Nh 

[10l9cm-9] 

^^total,Hotgas 

[1019 Mol 

FoVl 

179.83 

65.03 

2.59 

4.57 

FoV2 

17.73 

-52.25 

1.31 

2.31 

FoV3 

91.49 

47.95 

2.96 

5.21 

FoV4 

35.97 

-29.86 

2.81 

4.96 

FoV5 

289.95 

64.36 

2.62 

4.62 

F 0 V 6 

92.14 

-25.34 

2.06 

3.63 

FoV7 

287.46 

22.95 

2.62 

4.62 

F 0 V 8 

40.27 

-34.94 

2.59 

4.57 


TABLE 3 

Mock observations of eight randomly distributed 

DIRECTIONS IN THE SKY, IN OUR G.321 SIMULATION AT Z = 0. ThE 
total HOT GAS MASS IN THE SIMULATION, UP TO R„ir, IS 

M^ir = 1.2x10^° M0. We have defined hot gas as gas at 
T > 3x10® K. We have assumed L = 230.1 kpc = Rvir- 



EM 

[ 10 -®cm-® pc] 

L 

[kpc] 

n 

[ 10 -'‘cm-®] 

^^total,Hotgas 

[ 10 “ Mol 

FoVl 

2.03 

40.2 

2.2 

0.11 

FoV2 

0.16 

111.3 

0.4 

0.40 

FoV3 

1.58 

53.7 

1.7 

0.21 

FoV4 

0.79 

56.6 

1.2 

0.17 

FoV5 

0.89 

80.5 

1.1 

0.42 

F 0 V 6 

1.40 

59.1 

1.5 

0.25 

FoV7 

0.98 

74.1 

1.1 

0.36 

F 0 V 8 

0.61 

116.5 

0.7 

0.88 


TABLE 4 

Emission Measure (EM) values from the same eight 
RANDOMLY DISTRIBUTED DIRECTIONS IN THE SKY AS IN TaB. [S] 
Using EM and Nh values we have computed the path 
LENGTH (L) AND HOT GAS VOLUME DENSITY (n). FINALLY WE 
SHOW THE TOTAL HOT GAS MASS OBTAINED WHEN ASSUMING A 
SPHERE WITH CONSTANT DENSITY (n) AND ITS RADIUS EQUAL TO 
THE PATH LENGTH (L). ThE REAL TOTAL HOT GAS MASS IN THE 
SIMULATION IS M^ir = 1.2x10^° Mq. 


Planck Collaboration et al.| 2014) but are comparable 
to reported values tor isolated spiral galax ies based 


on weak leasing and stellar mass are 0.056 (Hoekstra 
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I [deg] 


Fig. 13.— Hot gas (T > 3x10® K) column density (top) in a full 
sky view of simulation G.321 in galactic coordinates, at z = 0. All 
values have been computed as observed from R = 8 kpc and at an 
arbitrary azimuthal angle. 


et al.|[2005). Such mismatch between cosmic and galaxy 


baryon budget suggests that either gas has escaped 
halos or never reached them and it and it is located 
in the Large Scale Structure warm -hot intergalactic 
(WHIM) gas (Cen & Ostriker 1999) or that warm-hot 
gas embedded in dark matter halos may account at 
least for part of the missing baryons if not all. In our 
simulations the baryonic fraction reaches (0.107—0.120) 
even after adding the hot gas mass component inside 
halos. Because we do not include the feedback effect of 
an AGN, that may eject the gas outside the halo, we can 
consider our estimates as an upper limit to the halo hot 
gas mass. Based on the previous discussion we argue 
that the missing baryons must be pl aced both like hot 


gas mass inside dark matte r halos (e.g. Gupta et al. 2012 


Miller & Bregman 

2013b) and also in the IGM, along 

the hlaments and 

ar from the main galactic systems 

(Rasheed et al.||2010 

Eckert et al.||2015). This is consis- 

tent with recent stu< 

res tor GGliVl in M31 dLehner et ai. 


5. CONCLUSIONS 

In this paper we have introduced a new set of MW- 
sized cosmological N-body plus hydrodynamics simula¬ 
tions. Their resolution and realism ar e high enough 
to allow us to st udy stellar kinematics (|Roca-Fabrega 
et al.||2013l |2014|) and evolution of large scale structures 
of these kind of galaxies. Most stellar disk parameters 
in our G.3 models are in agreement with observational 
ranges for the MW. Disk scale length and scale height 
as well as their dependence on the stellar age are re¬ 
produced. The simulated disk shows a large variety of 
large scale structures such us spiral arms, a ring and a 
bar. It also shows a flared and a warped configuration. 
The quality of our runs is comparable to th e mo st re¬ 
cent works prese nted by Guedes et al. (2011) and Mol- 
litor ei^al. (2015), a detailed comparison can be seen m 
Tablef^ Our simulations have a similar or a larger num¬ 
ber ofresolution elements (stellar particles or gas cells), 
higher spatial resolution, and a smaller mass per particle 
than the ones in those works. In addition, as we have 
shown in Fig.|^ our circular velocity curves are in a very 
good agreement with recent observations. On the other 
hand, a spheroid (bulge-|-halo) that is 2—3 times more 
massive than the one observed in the real MW is present 
in our G.3 models. Since the density of this spheroidal 
component is relatively low in the central region of our 
galactic systems it does not drive to centrally peaked cir- 
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Fig. 14.— Hot gas (T > 3x10® K) view in spherical galactic 
coordinates of the C.321 model, at z = 0. Top: metallicity as 
function of the azimuthal angle. Bottom: metallicity as function 
of the vertical angle. Zla and ZII are metals from SNe la and H, 
respectively. All values have been computed as observed from the 
galactic center. 

cular velocity curves. 

We have also analyzed the SFH in our simulations. We 
have found that our SFR values at z=0 are lower than 
recent observational values. Also the SFH peak in our 
G.3 model s occurs at higher re dshift than the one pre- 


dieted by (Behroozi et al. 

2013 

). In the most recent and 

realistic iVlW-sized Simula 

tions (Guedes et al. 12011 

Mol- 

litor et al.||2015 

Bird et al. 

2013), the predicted Sl'H is 


also not well reproduced. We have also found that total 
gas mass and its distribution in temperature slightly dif¬ 
fer from MW observations. This mismatch in SFH and 
mass gas distribution is usual in simulations due to the 
lack of understanding of subgrid physics and to the diffi¬ 
culty of accounting for all physical processes involved in 
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Fig. 15.— Hot gas (T > 3x10^ K) projected metallicity (color 
map) and galactocentric velocity direction (arrows) in the principal 
plane X-Y of the G.321 model, at z=0. Zla and ZII are metals from 
SNe la and II, respectively 
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Fig. 16.— Hot gas (T > 3x10® K) projected metallicity (color 
map) and galactocentric velocity direction (arrows) in two principal 
planes, X-Z (top) and Y-Z (bottom) of the G.321 model, at z = 0. 
Zla and ZII are metals from SNe la and II, respectively 


the formation and evolution of galaxies. Studies like ours 
will help to constrain such subgrid models. Recently, it 
has been proposed that the inclusion of a kinematic feed¬ 
back might solve the early star formation issue and thus 
change the position of the SFH peak. 

The work presented here is focused in the study of the 
distribution and the amount of potentially X-ray emit¬ 
ting gas (T > 3x10® K) of our simulations. The com¬ 
parison of the properties of the hot gas in our simulation 
with the observations of the hot halo X-ray corona of the 
MW and external galaxies was made by producing hot 
gas column density mock observations. Conclusions on 
this part follows: 

• Hot gas is not distributed homogeneously around 
the galactic halo. 

• There is no evidence of a hot gas thick disk in our 
models. 


• By using the real mass of warm-hot gas in our sim¬ 
ulations we conclude that about 53%—80% of the 
missing baryons can be accounted by hot gas in 
the halo corona. The exact value depends on the 
simulated galactic system (e.g. halo mass). We 
argue that the lacking baryons are placed in fila¬ 
ments (IGM), not inside main galactic system s, as 
preheated warm-hot gas (Rasheed et al. 2010). It is 
important to mention that in our simulations AGN 
feedback from a central galactic black hole has not 
been implemented. It is still open to debate how 
a relatively low mass black hole, as the one in the 
real MW, could affect the mass and temperature 
distribution of the hot halo gas and thus the SFH. 


• Hot gas density follows several streams in the halo 
that can be detected because of their changes in 
kinematics and metallicity. We have seen how low 
metallicity streams come from the IGM while oth¬ 
ers with higher metallicity depart from the disk. 
We have also detected the infall of high metallic¬ 
ity gas that can be interpreted as coming from a 
satellite accretion. 


• The total hot gas mass in the halo can not be re¬ 
covered when assuming a spherical uniform gas dis¬ 
tribution. In our model the real hot gas density 
distribution lead to an overestimation when using 
the mean column density plus a fixed optical depth 
under the spherical uniform assumption. On the 
other hand it drives to an underestimation of the 
total hot gas mass when we use the column density 
plus the emission measure when using the same as¬ 
sumption. 

• Following previous works we have found that a 
clear relation exist between total hot gas mass and 
virial mass of halos. This result is important as it 
becomes new method to constrain the total mass 
of galactic systems, such as the Milky Way. 
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Fig. 17.— Warm-hot gas slice plot at an arbitrary plane of the simulation G.323, at z = 0. 
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Fig. 18.— Circular velocity curves of G.321, G.322, G.323 and G.321_lr computed using GM/r approximation (see Tabj^for more details 
about models). Black: total, red: dark matter, blue: stars, green: gas. 
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APPENDIX 


A. RESOLUTION STUDY 

In this appendix we show a test of numerical convergence using a low resolution resimulation of our main model. 
The resimulation presented here (G.321_lr) has been obtained using the same IC and initial parameters as in G.321 
but with a resolution of 218 pc instead of 109 and ~9xl0® DM particles inside Rt, instead of ^7x10®. In Figfisl we 
show a comparison of the total, DM, stellar and gas circular velocity curves computed using y^GM/r) as a proxy 
for Vc- We see that resolution does not affect the total circular velocity curve profiles. However, it is interesting to 
note that circular velocity curve from the G.321 gas co mponent shows significan tly higher Vc values at large radii. 
This result is in agreement with what was obtained by Scannapieco et al. (20121; i.e. the properties of the gaseous 
component seem to be the most sensitive to numerical resolution ettects. We show in Tabj^ the values of the main 
parameters of the G.321 run and compare them with the corresponding ones of the G.321dr model. We see that, aside 
from the mass of cold gas and related q uantities, like SFR at z=0 , most parameters do not differ from each other 


by ore than 30%. This was also found in Scannapieco et al. (|2012|). From this result, we conclude that although, in 
general, parameters are not sensitive to changes in resolutiori, caution must be exercised when analyzing the cold gas 
component and the SFR. In addition, we notice that the cold gas fraction is not only sensitive to numerical resolution 
(the maximum refinement level) but also to changes in terms of number of cells (different aggressive refinement, see 
Sec.3.3). On the other hand, the difference observed in G.321 gas velocity curve at large radii is simply due to the fact 
that this run ends up wi t h much more gas. Additional stu dies of numerical convergence of the code can be found in 
Avila-Reese et al. (20111; Gonzalez-Samaniego et al. (2014). 
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G.321 

G.321 lr 

1 (G.321-G.321 lr)/G.321| X100 

Mvir [Mq] 

7.33x10“ 

7.33x1011 

0.0 

Tyir [kpc] 

230.1 

230.5 

0.2 

M* [Mq] 

e.ixioio 

6 . 0 x 101 ° 

2.0 

Mgas [M©] 

2.70x101° 

1.93x101° 

28.5 

^hotgas [M©] 

1 . 22 x 101 ° 

1.06x101° 

13.1 

"^warmgas [^o] 

5.66x10° 

7.33x10° 

29.5 

^coldgas [M©] 

9.34x10® 

1.37x10° 

85.3 

M 200 [M©] 

6.84x1011 

6.61x1011 

3.4 

1200 [kpc] 

175.6 

175.6 

0.0 

Fi),U 

0.120 

0.108 

10.0 

c 

28.5 

25.4 

10.9 

Rd [kpc] 

2.56 (4.89/2.21) 

2.13 (4.2/1.8) 

16.8 

^z,young [pc] 

277 

390 

40.7 

hz.oid [pc] 

1356 

1032 

23.9 

O'X 

-0.62 

-0.60 

3.2 

SFR (z=0) [M©yr-i] 

0.27 

0.16 

40.7 

Vc 0 (R= 8 kpc) [kms“^] 

239.8 

239.2 

0.3 

^peak [kpc] 

5.69 

5.69 

0.0 

VclRpeafc) [kms-l] 

243.8 

243.6 

0.1 

V 2 . 2 /V 2 OO 

1.90 

1.87 

1.6 


TABLE 5 

Comparison between general parameters of GARROTXA original model G321 and low resolution run G321_lr. All 

PARAMETERS SHOWN HERE HAVE BEEN WELL DESCRIBED IN TAB^ 


B. JEANS LENGTH 


In this section, following the condition proposed by Truelove et al. (1997), we show that in our simulations the 
resolution is above the minimum that ensures the formation of stellar particles is of physical origin rather than 
numerical. In order to see this we have computed the local Jeans wavelength of the cold gas (Aj) according to Eqj^ 
In Figj^we show a histogram of the Jeans length in units of the size of the corresponding cell for simulation G322. 
Except lor few ones, cells satisfy the Truelove criterium. 
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Fig. 19. 
simulation 



Aj/Cell size 


— Histogram of the cold gas (T<4xl0'^) Jeans length computed following EqlH in units of the corresponding cell size, for 
G322. 





